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I Al)^^ract 

I We  have  studied  the  Brillouin  scattcriii};  spectra  of  oligomers  of 

polyethylene  glycol,  along  with  ethylene  glycol.  These  systems  vary 
greatly  in  viscosity  as  well  as  molecular  weight.  We  have  found  that 
except  for  ethylene  glycol  (the  monomer),  the  Brillouin  frequencies 
and  the  linewidths  of  the  rest  of  oligomers  show  more  or  less  the  same 
temperature  dependence,  independent  of  viscosity  and  molecular  weight. 
These  results  indicate  that  the  motion  responsible  for  the  hypersonic 
dispersion  and  attenuation  is  localized.  We  have  also  developed  a 
theory  based  upon  the  Zwanzig-Mori  formalism  to  account  for  the  ex- 
perimental results.  We  have  then  shown  that  in  order  to  obtain  a 
good  agreement  between  theory  and  the  experiment,  the  usual  Markov 
approximation  describing  the  dynamic  behavior  of  the  memory  function 
cannot  be  used.  Thus,  the  usual  hydrodynamic  ec\uations  with  frequency 
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1 lit  roiliict  ion 

Modern  development  in  optical  technology  has  made  the  Rayleigh 
and  Brillouin  light  scattering  technique  very  useful  for  the  investi- 
gation of  molecular  dynamic  properties  of  large  molecules  sucli  as  poly- 
mers.^*"* In  general,  the  Rayleigh  and  depolarized  Rayleigh  measure- 
ments yield  information  on  the  translational  and  rotational  diffusive 
motions  of  molecular  segments  and  an  entire  molecule.'  On  the  other 
hand,  Brillouin  scattering  is  due  to  tlie  interaction  of  incident 
light  with  thermally  induced  sound  waves  in  the  medium  and  thus  re- 
flects the  cooperative  motion  among  different  segments  and  molecules.^ 
Therefore  the  interaction  between  different  molecules  and  segments  is 
expected  to  have  an  important  effect  on  the  Brillouin  scattering 
spectrum  of  a polymer  fluid. 

We  have  previously  reported  the  Brillouin  scattering  results  of 
polyriropylene  glycol  Prom  these  experimental  results,  we 

have  come  to  believe  that  the  molecular  dynamics  which  can  affect 
Brillouin  spectra  of  polymeric  liquids  are  localized  motion  involving 
only  very  few  segments.  In  order  to  further  affirm  these  results,  we 
have  carried  out  a systematic  investigation  of  the  Brillouin  spectra 
of  low  molecular  weight  oligomeric  polyethylene  glycol  (PnC)  liquids 
including  ethylene  glycol.  The  experimental  setup  for  such  studies 
was  similar  to  that  used  for  the  PPH  measurement.^’^  We  have  also 
carried  out  viscosity  measurements  and  compared  the  Brillouin  scat- 
tering results  with  the  classical  theory  on  the  sound  propagation  and 
attenuation.  We  have  found  tliat  classical  theory  is  not  adeejuate  for 
PKCi.  The  viscosities  of  bulk  polymer  fluids  such  as  PP(i  and  Pl-G  are 
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much  higher  than  those  of  ordinary  sim|)Ie  li(|uids.  This  is  due  to 
I strong  intermolecular  interactions  among  segments,  which  render  these 

polymer  fluids  highly  viscous.  In  a previous  paper^  we  have  used  the 
Zwanzig-Mori  linear  response  theory  to  analyze  the  Brillouin  scatter- 
ing results  of  bulk  liquid  polypropylene  glycols  (PPfi) . It  was 
shown  that  the  second  moment  of  the  intersegmental  interaction  poten- 
tial and  its  relaxation  dynamics  have  a dominant  effect  on  the  line- 
width  and  fi'cqucncy  sliift  of  the  Brillouin  peak. 

In  the  present  paper,  we  carry  out  a similar  analysis  for  poly- 
ethylene glycol  (PI-G)  liquids.  To  our  knowledge  the  experimental 
Brillouin  scattering  results  of  these  fluids  have  not  been  reported 
in  the  literature.  We  have  presented  the  theoretical  result  more  com- 
pletely here  than  in  the  earlier  paper.  The  derivation  is  carried  out 
in  terms  of  two  coupled  variables  (density  and  velocity;  or  their  lin- 
ear combinations).  Under  cetain  approximations  the  final  result  in  the 
present  paper  reduces  to  the  previous  one  (Hq . (26)  of  Ref.  5),  with 
which  the  PPG  Brillouin  scattering  results  were  analyzed.  The  simpli- 
fied equation  is  used  to  analyze  the  experimental  Brillouin  scattering 
results  of  PEG.  The  least  souares  fitting  of  PEG  NIIV  200  to  the  theo- 
retical expression  was  carried  out.  The  result  was  found  to  be  simi- 
lar to  tliat  of  PPG  MW  425  as  shown  in  tlu'  pvevicais  paper. 

'llic  mc.asured  sound  velocity  .-ind  attenuat  i mi  coefficient  results 
of  the  Brillouin  scattering  spectra  of  PEG  dimer,  trimer  and  oligomers  i 

of  average  molecular  weights  of  200  and  400  are  shown  together  witli  | 

that  of  ethylene  glycol  monomer  in  Figs.  1 and  2.  The  results  indi-  ; 

cate  that  Brillouin  data  show  significant  dispersion  near  200  K. 

The  Brillouin  scattering  result  of  the  monomer  is  quite  ilifferent  from 
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tl\e  former  uroup  of  I’l.Ci,  wliose  Ui'illoiiin  scattering  results  arc  found 
virtually  independent  of  molecular  wciKltt.  Ihe  cause  for  tliis  differ- 
ence is  discussed.  We  have  compared  the  Brillouin  scattering  results  of 
the  PPfi  with  the  Pl;(l  series,  and  have  found  a significant  difference  be- 
tween them.  This  difference  is  interpretcil  as  due  to  their  different 
segmental  mass  and  different  degree  of  i ntermolecular  interactions. 

The^rt^l  1 Cons  i de r.-U  i on 

For  a polymer  liquid  made  up  of  a collection  of  identical  chain 
segments,  the  isotrojiic  part  of  the  polarized  Bri 1 louin-Ray lei gh  spec- 
trum observed  at  scattering  angle  6 is  determined  by  the  Fourier  trans- 
form of  the  time  correlation  function  of  density  fluctuations,  given 
by 

Cpp(<l.f)  = <6pfq,t)  fio  (q)>  (1) 

where  6p(q,t')  is  the  qth  mode  of  the  fluctuation  of  density  and  is 
given  by 

6p(q,t)  = T.  (2) 

j 

Here  r.(t)  is  the  position  of  scatterer  (\^e  take  it  to  be  a chain  seg- 
ment) at  time  t,  q is  the  scattering  wave  vector  and  is  related  to  the 
wave  vectors  of  incident  and  scattering  light  by 

q = Ks  - ^i  (.^) 

In  F(i.  (2)  the  summation  over  j is  to  include  all  scatterers  within 
the  scattering  volume. 

In  order  to  describe  the  Brillouin  spectrum,  we  need  to  consider 
V(q,t)  in  addition  to  the  density  fluctuation  6p(q,t).  Here  V(q,t)  is 


given  by 


Vf(i,t)  = V.  (4) 

.i  '■ 

where  Vj  is  the  velocity  of  segment  j.  I f we  choose  q to  be  in  the 
direction  of  the  z-axis  of  the  laboratory  coordinate  system;  then  by 
symmetry  only  the  z-comjionent  of  VCq.f)  can  couple  to  6p(q,t).  In 
principle  we  need  to  consider  the  energy  fluctuation  to  account  for 
the  dynamics  of  density  fluctuations  jiroperly,  but  tlie  inclusion  of 
the  energy  fluctuation  can  only  affect  the  Brillouin  spectrum  in  minor 
ways  (one  of  which  is  to  render  tl»e  liypersonic  wave  to  propagate  with 
an  adiabatic  velocity,  rather  than  with  an  isothermal  velocity,  as  con- 
sidered in  this  paper).  A comiilote  microsconic  theory  on  Brillouin 
scattering  including  the  energy  fluctuation  is  presently  in  progress. 

The  qth  mode  of  density  fluctuations  when  coupled  to  V^(q)  will 
give  rise  to  two  Brillouin  peaks  in  the  Br i 1 louin-Rayleigh  spectrum. 

The  two  are  the  Stokes  and  anti -Stokes  peaks  associated  with  the  lon- 
gitudinal hypersonic  wave.  As  shown  in  a previous  paper,  the  two 
Brillouin  peaks  can  be  calculated  by  using  the  Zwanzig-Mori  equation^ 
using  two  dynamic  variables  and  C , which  are  defined  as 

where  for  the  simplicity  of  notation  we  have  s.uppressed  the  arguments 
associated  with  ^ , Ap  and  , understanding  that  thev  all  dcjiend  on 
both  (]  and  t. 

It  is  easy  to  show  by  using  time  reversal  symmetry  that  at  t=o , 

* ♦ * * 

f,  > = 1 and  > = ''f.  > = o. 


In  the  previous  paper,  we  nej^lected  the  dynamic  couplinj^  between 


> anil  and  treated  only  the  motion  of  and  f as  inde/)endent  vari- 

ables. In  this  paper  we  take  this  coupling  into  account,  and  for  this 
purpose,  we  define  a column  vector  A consisting  of  and  f,  as  its 
elements 


The  equation  of  motion  for  Aft)  is  given  formally  by  the  Liouville 
equation , 

Aft)  = iLA(t)  (7) 

where  L is  the  classical  Liouville  operator  governing  the  time  evolu- 
tion of  the  dynamic  vector  variable  Aft).  While  it  is  straightforward 
to  write  down  a formal  e.xpression  for  I,  in  terms  of  the  kinetic  and 
potential  energy  operators  in  a many  body  system,  in  practice  this  is 
seldom  done  in  t polymer  system  due  to  its  large  number  of  internal 
degrees  of  freedom  and  the  intermolecular  and  intramolecular  potential 
functions.  Thus  we  shall  utilize  only  the  formalistic  symmetry  prop- 
erties of  the  Liouville  operator  such  as  the  time  reversal,  reflection 
and  translation  symmetry  properties,  instead  of  writing  out  its  ex- 
plicit form.  In  its  formal  definition  the  Liouville  ojierator  is  self- 
adjoint  . 

F'or  the  set  of  dynamic  variables  Aft),  the  •'tori  transformation 
of  l-q.  f7)  gives. 

^ Aft)  = i^Aft)  - dr  KfT)-Aft-T)  + Fft) 

(8) 


i 


I 


where  the  frequency  matrix  S)  is  );iven  in  terms  of  I. 


Q = <(1.A)  A > 


The  otlier  (]uantitics  are 

Kfx)  = <I'(t)  f''> 


FCt)  = e 


iQl.T, 


(9^ 

flO) 

(in 


F = iQI.A  fl2) 

where  Q = (l-p),  )i  hoini;  the  Zwanz i p-Mor i projection  onerator  de- 
fined as 

pt;  = <(;a'^>  <aa''>'1a  = <('.a^>a  (13) 

due  to  the  fact  that  <AA^>  is  a unit  matrix. 

Substituting  liq . (fi)  into  F.t| . (Oj  , we  obtain  in  the  limit  of 
sma II  q , 


SI  = 


( ID  O 

O 


O -ID 


(14) 


where 


?r  2 


= q"^  q^(Dmx,p)  ‘ = q'^C.j.^.  Here  Xy  is  the  iso- 


thermal compressibility  and  is  related  to  <|6()|^>  by 
<|iSp|^>  = Vp2kTx.p 

as  q — >■  0.  is  the  isothermal  sound  velocity.  l 


(15) 


liquation  (14)  thus  indicates  that  A(t)  jiropaRates  with  the  iso- 
thermal sound  velocity,  instead  of  the  desired  adiabatic  sound  veloc- 
ity. The  reason  for  this  is  that,  as  mentioned  .above,  we  h.ive  not  in- 
cluded the  enerK)'  fluctuation  term  in  the  calculation.  Since  we  sluU  1 
at  the  end  treat  as  a parameter  to  fit  the  experimental  data,  the 
omission  of  the  coupling  of  density  to  energy  will  not  affect  the 
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interpretation  of  the  Brilloiiin  ilata  to  he  niven  in  the  next  section. 

The  matrix  K(t)  is  known  as  the  memory  function  matrix.  Upon 
the  substitution  of  hqs.  (5)  and  ((>)  into  Ii(|.  (10)  and  after  carry- 
ing out  tlie  necessary  matrix  multiplication,  we  find  for  the  memory 
function, 

^K(t)  -K(t) 

-K(t)  K(t) 


K(t) 


(16) 


where  for  small  q 


K(t)  = jf<|e*'‘^(il,V)l  fil.vf>  - u "1 
/ () 


(17) 


i I.T  iQl.T 


For  brevity  we  have  used  V for  V , and  c ' for  e 

z 


The  latter 

simplification  is  possible  for  q — ► ().'  At  i = 0 it  has  been  shown 
in  Appendix  II  of  Ref.  .S  that  K(0)  in  the  small  (]  limit  has  the  ana- 
lytical form: 


=1  ^‘2Ni  ^ ^ <vr-r--.~w  0=’-  - Ui  2} 


.-^kT 


3(Z.-Z.)2  '1  -J 
1 .1  ^ 1 J 


(r-  - 0) 
m q 


(18) 


where  U is  the  total  potential  energy  of  the  system.  The  first  term 
of  Eq.  (18)  can  be  shown  to  be  positive  and  greater  than 
can  easily  recognize  that  this  term  is  proportional  to  the  spatial 
second  moment  of  the  intermolecular  and  intersegmental  ]iotential 
energy  in  the  Z direction. 

Substitution  of  Eqs.  (6),  (14),  and  (16)  into  the  generalized 
I.angevin  e(|uation  (luj.  (8))  gives  two  coujiled  kinetic  ecpiations: 
ar,  I 

■g:^  (t)  = ± (i)^f,^(t)  - fK(T)f,^  (t-T)  - K(T)f,  j't-T)ldT  + F^(t) 

(ID) 
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where  are  the  fluctuation  forces  actinjj  on  tlie  slow  variables 

f ^ ( 1 1 . I'he  two  coupled  kinetic  ecjuations  can  be  solved  by  the  method 
of  l.aplace  transforms. 

The  memory  function  K(t)  is  tl\e  tim(‘  correlation  function  of 
fluctuation  forces  and  sliould  always  decay  to  zero  as  time  goes  on. 

We  have  discusseil  in  Uef.  .S  that  the  time  depciulence  of  K(t)  can  be 
approximated  as  an  exponential  decaying  function  with  correlation 
time  Tjj, 

Kfr!  = KfO)  e (20) 

we  denote  f.^(s)  and  K(s)  as  the  b.aplace  transforms  of  f^(t)  and  K(tJ 
respectively,  and  also  denotes 

I,(s)  = <(:,(s)  c/(0)>  (21) 

( fsl  = <f  fs)  *(0)>  (22) 

From  l;<i.  (1)  and  the  definitions  of  (Ihl.  (5))  one  can  show  that 
the  light  scattering  spectrum  due  to  density  fluctuation  is 

1(a))  = Re  [I^fs)  + I_(s)  + I^_(s)  + I_^(s)l 

s— > io) 

where  Re  denotes  taking  the  real  part  after  s is  substituted  by 

s — >•  io) 
io). 

.Solution  of  Fi],  (19)  using  FRi.  (20)  for  K())  yields  the  following 

1(a))  = T^fo))  + I (o))  + (o))  + I ^(o)) 

M - A, 

• ) = fM  v'aVi'^  V [,r*ar~r N "--iT'p 


(2.S) 


(24) 


results : 
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Y,  = f {N  - (31) 

where  M ;irnl  N are  the  rial  and  ima^jiiiary  parts  of  Kficii),  respec- 
tively. In  reaching  the  results  of  lai.  fJ.S),  we  have  neglected  an  un- 
important term  which  appears  in  the  denominators.  Since  the  1 Ci.i) 

■*  T 

terms  contribute  much  less  to  the  total  spectrum  than  the  l,(<ii)  terms, 
the  neglect  of  the  unimportant  term  in  1 ^ (h\)  is  justified  as  far  as  the 
calculation  of  the  total  Brillouiti  spectrum  is  concerned. 

In  the  temperature  region  where  KT|^^<<1  is  valid,  and  torms 
of  F.q.  (24)  are  much  smaller  than  the  other  terms,  and  I^(a))  gives  a 
maximum  spectral  intensity  around  m = + |n|).  From  Eq . (25) 

one  notices  that  1^  (oi)  is  identical  to  1 ^(m).  They  give  rise  to 
a spectrum  with  a much  broader  bandwidth  than  that  from  T^(ut).  These 
terms  have  an  intensity  maximum  centered  at  w = 0.  However,  the  spec- 
tral wings  of  (w)  in  the  v i c ini  ty  of  w = ± ( | 1 + |n|)  are  negligi- 
ble compared  with  l^((ii). 


The  1^  ((ji)  terms  as  given  in  licj.  (25)  rejiresent  the  dyn.'imic 
coupling  of  tf)e  two  phonon  inoiles  with  opposite  velocity.  This  term 
is  important  to  account  for  the  shane  of  the  entire  Br i 1 louin-Ray- 
leigt)  spectrum,  as  it  gives  rise  to  a spectrum  I'entercil  at  zero  fre- 
quenc>',  as  [lointed  iiiit  in  t lu'  ('arli(>r  paper.'’  However,  despite  the  pi 
sence  of  fi.O  , we  have  found  th.it  I (m),  with  tlie  A and  H teniis  nt 

t • ♦ ♦ 

glected,  jtioviilos  an  adc()uatc  account  for  tlu'  temperature  dependences 
of  tite  Brillouin  linewidth  and  tlie  frequency  shift.  The  simplified 
spectral  functions  as  given  by  the  1 ^ (wl  terms  were  used  to  analyze 
the  Brillouin  data  of  PPG  in  the  earlier  paper.  We  now  use  the  same 
spectral  function,  coupled  with  the  least  sipuircs  fitting  [irocedure 
to  aisilyze  the  I’l.G  Brillouin  sc.ii  tei'i 'm>  results. 


Discussion  of  the  iixper imenta  1 


llesults  of  Pr.Ci  J^riJJouin  .Scattering 


The  ne.isured  Brillouin  shift  ( in  Hertz)  is  related  to  the  hy- 


personic velocity  by  the  equation: 


V = f A /_  . ,0, 

s Bo  2nsin(Ty) 


where  is  the  wavelength  of  the  incident  light  in  vacuum,  n is  the 

index  of  refraction  of  the  scattering  medium.  0 is  the  scattering 
angle  from  the  incident  light.  The  attenuation  coefficient  over 

one  wavelength  of  the  hypersound  is  expressed  as 


where  is  the  instrument  corrected  Brillouin  spectral  full  width  at 
h.'il  f-max  imum  intensity.  is  the  wavelength  of  the  hypersonic  wave 
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and  is  {{iven  by 


A 

s 


f34) 


1 lu'  expo  I- i men  t a I results  of  .iiul  of  eth)'lene  glycol  (mono- 
mer) and  I’iU:  dimer,  t rimer,  ami  i>oIymers  of  MK  2(10  and  Ml\'  400  at  Do” 
scattering  angle  are  presented  in  F-igs.  1 and  2 as  functions  of  temper- 
ature. fxcept  for  the  ethylene  glvcol  results,  the  nX  and  V results 

s s 

^or  t’'e  otl'or  olinopcrs  are  virtuallv'  indenondent  o*^  molecular  weight. 
As  nointed  oat  in  t he  i ntro'’uct inn . the  sane  molecular  weieht  indcoen- 
dence  was  e'so  observed  “^or  tiu'  I’l’b  series.^  Since  the  Brillouin  scat- 
tering spectrum  of  a viscous  fluid  is  greatly  influenced  by  the  inter- 

molecular  and  interscgmontal  interactions,  the  distinctive  aX  results 

s 

of  ethylene  glycol  indicate  that  the  dynamics  of  the  intermolecular 
interaction  of  the  monomer  is  quite  different  from  the  other  members 
of  the  series.  The  difference  is  most  likely  due  to  its  short  molec- 
ular dimension  and  the  presences  of  dominant  hydrogen-bonding  inter- 
action with  neighboring  molecules. 

Despite  their  different  molecular  weights,  we  have  found  that  all 
of  the  I’HG  samples  under  the  jiresent  study,  have  practically  the  same 
index  of  refraction  (to  within  ]%  of  error).  Thus,  according  to  F.qs. 
(.■^2)  and  (.A.’S)  , the  same  index  of  refraction  implies  that  the  Brillouin 
frequency  shift  and  the  linowidth  are  practically  the  same  for  all  of 
t lu'  I’lCi  .imples  of  diffh'rcnt  moh'cular  weight'-. 

According,  to  t lu'  Na  v j er -S  t okes  ('i|u.i  t i on  , the  hypersonic  velocity 
depends  only  on  tidiabatic  compressibility,  and  the  attenu.ition  coeffi- 
cient (or  linewidth)  otily  on  the  static  shear  and  hulk  viscosities,  for 


the  I’lid  system,  the  static  shear  viscosity  clianjjes  markedly  with  tem- 
perature; the  static  hulk  viscosity  which  is  more  or  less  proportional 
to  the  shear  viscosity  also  clianj,;cs  rapidly  with  temperature.  Thus 
the  classical  theory  of  sound  projiajjat  ion  and  attenuation  cannot  he 
used  to  describe  the  liypersonic  propa(>ation  and  attenuation  behavior. 

It  is  now  recoftnized  tliat  tlie  .acoustical  beliavior  of  highly  vis- 
cous li(|uiils  ili-|)iiid:.  on  tiu'  product  of  snimd  fretpiency  and  the  viscos- 
ity (i.e.,  Ij^n)  and  not  on  the  p.irameters  sop,uatel\.  lor  tliis  reason, 
we  have  me.isured  the  viscosity  of  ethylene  glycol  and  all  of  the  mem- 
bers of  I’liC  studied.  We  have  jilotteil  the  hy]H’rsonic  velocity  (V)  and 
attenuation  coefficient  (aA)  as  a function  of  fj^n  in  Figs.  (3)  and  (4) 
respectively.  In  terms  of  such  plots,  V and  itX  separate  into  differ- 
ent curves  for  different  I’bC  members,  due  to  the  fact  that  each  mem- 
ber of  ri;G  has  a different  viscosity  ,at  a given  t cmper.ature . A.t  a 

fixed  f^n,  the  monomer  gives  the  largest  V.  However,  the  fiinction.al 
n 

dependence  of  nX  on  f„n  is  more  complicated.  At  low  f_n  value  ( fr,n 
< ISO  GHz  cp) , corresponding  to  the  high  temperature  situation,  ntX 
is  independent  of  molecular  weight,  beyond  this  value  c«X  i net  oases 
with  increasing  degree  of  polymerization,  witli  some  anomaly  observed 
for  the  monomer  result.  The  velocity  data  approach  plateau  values 
at  large  fj^n  for  each  oligomer,  whereas  nX  decrc.ises  gradually  at 
largo  fj^n. 

Houbtless  to  say,  the  f|^n  dependence  of  V and  oX  is  due  to  the 
structural  relax.ation  effect  taking  place  at  the  hypersonic  fretpiency 
region  in  I’fG.  Assuming  tluit  the  structur'.'il  relaxation  time  is  propor- 
tion.il  to  till'  shear  viscosity,  we  ha\e  tried  to  interpret  our  ri'sults 
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ji 

{\ 
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iKH'oi'd  i nj;  to  a rolaxatioii  llu'oiy  set  t'orlli  by  I ('out  ov  i ill , ^ aiul 
I ouikI  that  till'  relaxation  tlieoi)'  |)tii.liits  too  stroiij’,  a ilcMiciuli'iKC 
of  V aiul  oA  on  fj^ri  at  large  f^'i. 

We  thus  fit,  with  the  help  of  a h'aat  squares  [irograin,  the  cx- 
[ler i men t a 1 results  of  the  frequency  shift  .iml  linewidth  to  the  values 
generated  from  the  spectral  functions  for  Hrillouin  scattering  ;is 
given  in  fq.  (24).  Neglecting  the  and  terras,  we  obtain  from 

f-q.  (24), 


I(ra)  = 


M? 


_ _M_ 

fw-w 


-NT?  ^ M2' 


M 

[ (ij+oj  -N  I? 


(.^.S) 


In  the  1 oast  - S(|ua  res  fitting,  we  first  g.enerated  the  whole  Brillouin 
spectrum  using  hc).  (.^5).  The  froipiency  shift  of  ])cak  maximum  and 
linewidth  of  the  generated  spectrum  were  then  read  and  compared  with 
the  experimental  values.  The  difference  between  the  theoretical  and 
experimental  values  is  then  minimized  with  the  1 east -s()uares  fitting 
program.  Ihe  temperature  depeiulence  of  the  relaxation  time  t|^  asso- 
ciated with  the  memory  function  is  assumed  to  follow  the  Arrhenius 
eipiat  i on 


where  is  the  activation  energy  for  the  relaxation  process  of  the 
second  moment  of  i ntersegmenta  1 interaction  potential.^  In  the  le.ast- 
stpiares  fitting  program,  there  are  three  adjustable  parametci's:  the 
initi.il  value  of  memory  function  .at  t = 0,  K;  ttie  activation  energy 
li.i  ; and  the  relaxation  time  1 ,at  infinitely  high  temiieratui'e  (see 
lii] . (.'^b)).  The  result  of  the  1 e.ast -sipi.ares  fitting  is  shown  in  fig. 
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5,  alDtij.;  with  t lu'  i men  t a 1 results. 

The  fittings  were  carried  niit  in  two  cases:  1)  independent  of 
temperature  and  2)  a function  of  temperature.  The  latter  case 
represents  a more  realistic  situation.  With  fixed,  the  fitting  to 
the  linewidth  data  is  excellent,  however,  there  are  deviations  in  the 
calculated  frequency  shifts  from  their  experimental  data  at  both  high 
and  low  temperatures.  As  pointed  out  in  the  previous  paper, ^ m , be- 
ing  proportional  to  (max.|.)  is  temperature  dependent  and  is  larger 
at  lower  temperature.  If  the  temperature  dependence  of  is  included, 
tlte  calculated  frequeticy  shift  values  are  significantly  improved  at 
both  low  and  high  temperatures.  Since  the  functional  form  of  the  tem- 
perature dependence  of  is  not  available,  to  demonstrate  this  point, 
we  simply  assume  that  the  temperature  ilependcnt  function  for  is 

11^  I in  2 m 

w =10  + — + — + — csyi 


q q 


T t’’ 


where  is  the  value  of  at  infinite  high  temperature,  and  ni]  , 
m? , and  m3  are  adjustable  parameters.  Using  f.q.  (.'^7)  , we  have  found 
that  the  fitting  to  the  linewidth  data  is  virtually  unchanged,  but 
the  fitting  to  the  frequency  shift  data  becomes  excellent  at  all  tem- 
peratures. The  improvement  of  the  fitting  to  the  frequency  by  allow- 
ing w to  be  temperature  dependent  does  not  change  the  K,  T.  and  x ” 

'I  tl  K 

values  significantly.  The  obtaim-d  host  values  for  these  three  para- 
meters are  K = 13  (GUz)^ , H.  = 4.11  kcal/mole  and  r,°  = 1.56  x lO'^'* 

n K 

sec  for  u fixeil,  and  K = 14.2  (Ullz)^;  f = 3.58  kcal/mole  and  i ° = 
q a K 

3.0  X 10‘''*  sec  for  being  temperature  dcnetident . The  difference  in 
the  values  of  activation  energy  is  l(":s  than  15b,  which  is  about  the 
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order  of  the  usual  error  in  the  activation  eiu‘r}',y  iiu'asurenu'nt  . 

As  in  the  least-squares  fitting  of  the  I'Pti  data,  the  obtained  K 
value  is  a fraction  of  • However,  the  activation  energy,  K , is 
found  slightly  smaller  than  that  of  PPG.^  Qualitatively  this  is  ex- 
pected, since  PPG  has  an  extra  methyl  group  for  the  side  chain,  which 
will  incrase  the  hindrance  against  the  motion  of  each  segment.  Never- 
theless, the  general  characteristics  of  the  results  obtained  for  PFiG 
by  the  least-squares  fitting  are  similar  to  those  of  PPG.  Thus  the 
discussion  of  tlie  previous  paper'’  on  PPG  applies  equally  well  to  tl)e 
tu'esenl  c f:<'. 

A common  feature  of  the  I’PG  series  studied  in  the  previous  paper 
and  the  PPG  scries  of  the  present  study  is  their  large  viscosities, 
compared  with  fluids  composed  of  simple  molecules.  In  Pig.  1 of  the 
previous  paper,  we  have  shown  that  the  hrillouin  scattering  results 
of  the  PPG  samples  are  virtually  independent  of  their  viscosities. 

This  is  also  the  case  for  the  I'PG  series.  The  Itigh  viscosity  of  the 
polymer  fluid  invalidates  the  Navicr-Stokes  equation  used  for  the 
Newtonian  fluids.*^  As  seen  from  our  theoretical  analysis,  the  memory 
function  is  characterized  by  the  spatial  second  moment  of  interaction 
potential  between  segments  and  its  relaxation  dynamics.  Hue  to  the 
strong  i ntersegmenta 1 interaction,  the  memory  function  cannot  lie  re- 
placed by  a relaxation  parameter  multiplied  by  a delta  function  gener- 
ally kiK'wn  ns  the  Markov  approximation.  Fn  the  present  tFieory,  the 
convolution  form  of  the  memory  function  and  tlie  dyn.imic  variables  is 
preserved.  I'his  feature,  although  not  well  neognized,  is  lulieved  to 
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be  very  important  for  describing  the  structural  relaxation  process 
j taking  place  in  the  hypersonic  freciucncy  region.  To  be  more  specific, 

the  imaginary  part  of  the  Fourier  transform  of  the  memory  function 
(the  N function  given  in  Fq.  (29))  is  very  important  to  account  for 
the  dispersion  of  the  peak  frequency,  and  the  real  part  (M  as  given 
in  hcj.  (28))  for  tlic  maximum  in  ttie  linewidtii. 

In  Frenkel's  hole  theory  of  liquids,'^  molecules  are  considered 
to  move  continuously  into  the  vacancies  which  are  proposed  to  exist 
in  a viscous  liquid.  Based  on  tliis  hole  concept,  Isakovicli  and 
Cliaban'*'  have  modeled  the  highly  viscous  liquid  as  a micro- inhomo- 
geneous medium  to  explain  the  sound  attenuation  at  tlie  hypersonic 
frequency  region.  A micro-inhomogeneous  medium  is  one  which  consists 
of  two  oi  more  inhomogeneous  regions,  the  dimensions  of  which  are  as- 
sumed to  l)e  small  compared  to  the  wavelength  of  sounil.  They  have  as- 
sumed that  upon  the  passage  of  sound  waves  an  exchange  of  energy  be- 
tween different  regions  in  the  micro-inhomogeneous  medium  occurs.  The 
energy  exchange  is  considered  to  he  due  to  mutual  diffusion  of  holes  in 
I different  regions.  Our  present  theoretical  analysis  does  not  assume 

I such  a model.  It  is  shown  that  the  initial  value  of  the  memory  func- 

I 

I tion  is  characterized  by  the  spatial  second  moment  of  the  intersegmen- 

tal  interaction  potential.  The  rapid  notion  of  a chain  segment  from 
one  position  to  another  is  responsible  for  the  decay  of  the  memory 
function,  and  also  causes  a local  reorganization  of  the  structure  in  a 
polymer  licpiid.  This  microscopic  i)icturc  is  implied  in  Isakovich  and 


Chahan's  model  of  stnicturni  relaxation  in  a highly  viscous  liquid, 
hut  they  do  not  include  the  memory  effect  in  the  calculation.  The  pre- 
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sent  theory  is  developed  on  tlie  basis  of  the  Zwanzig-Mor i formalism, 
and  allows  one  easily  to  relate  the  phettomenoloj^ical  quantities  in 
terms  of  a microscopic  picture  associated  with  the  physical  system 
under  study  with  a minimum  of  sulisiiliary  assumptions  about  the  na- 
ture of  the  system.  Tlie  present  microscopic  theory,  thus,  serves  to 
provide  the  Isakovich  and  Chaban  phenomenological  theory  with  a more 
satisfactory  statistical  foundation. 

Summary  and  (!onclus  ion 

We  have  studied  the  Brillouin  spectra  of  Pf.G,  along  with  ethylene 
glycol.  These  systems  along  with  PI’G  vary  greatly  in  viscosity  as  well 
as  molecular  weight.  We  have  found  that  except  for  monomer,  the  Bril- 
louin  frequency  and  the  linewidth  of  PPG  are  insensitive  to  viscosity 
and  molecular  weight.  These  results  indicate  that  the  motion  respon- 
sible for  the  hypersonic  dispersion  and  attenuation  is  localized.  De- 
spite some  difference  with  respect  to  the  temperature  behavior  in  the 
relaxation  time  associated  with  memory  function,  the  PPG  results 
closely  resemble  those  of  PPG. 

We  have  developed  a theory  on  the  basis  of  the  Zwanzig-Mori  formal- 
ism to  account  for  the  experimental  results.  We  have  then  shown  that 
in  order  to  obtain  a good  agreement  between  theory  and  the  experiment, 
the  usual  Markov  approximation  describing  the  dynamic  behavior  of  the 
memory  function  cannot  be  used.  Thus,  the  usual  hydrodynamic  equations 
with  frequency  independent  transport  coefficients  are  inadequate  in  de- 
scribing Brillouin  scattering  of  viscous  liquids.  Further  studies  of 
Brillouin  scattering  spector  will  definitely  provide  useful  informa- 
tion about  the  dynamic  structure  of  the  liquid  state. 
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Figure  1.  Hypersonic  velocities  of  ethylene  glycol  (o)  aod  PFG 
dimer  (A),  trimer  (<>1.  f'W  200  (#)  and  MIV  400  (□) 
as  a function  of  temperature. 

Figure  2.  Hypersonic  attenuation  coefficients  as  a function  of 
.temperature  (symbols  are  the  same  as  those  in  Fig.  1) . 

Figure  3.  Hypersonic  velocities  as  a function  of  f n (f,,  is  the 

B B 

hypersonic  frequency  and  n the  shear  viscosity;  other 
symbols  are  the  s;ime  as  those  in  IMg.  1). 


Figure  4.  Hypersonic  attenuation  coefficients  as  a function  of  f n 

B 

(symbols  are  tl\e  same  as  those  in  Fig.  3). 

Figure  5.  Comparison  of  the  calculated  and  the  observed  Brillouin 
peak  frequency  and  linewidth  data  as  a function  of  tem- 
perature. The  solid  line  is  obtained  for  being  tem- 
perature dependent;  the  daslied  line  for  fixed  w . 
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